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INTRODUCTION 


x y  Prccipi tation  static  generated  on  aircraft  by  corona  discharges 
has  long  been  recognized  as  a  major  problem  in  the  design  of  noise-free 
radio  communication  systems.  However,  little  attention  has  been  given  to 
the  radiated  electromagnetic  fields  that  accompany  the  discharge  process. 
This  report  describes  work  performed  by  the  Denver  Research  Institute  in 
the  investigation  of  the  electromagnetic  fields  radiated  by  conducting 
bodies  in  corona. 

Corona  discharges  occur  in  the  strong  local  electric  field  which 
builds  up  around  the  edges  of  an  airframe  as  the  airframe  potential  rises 

j 

from  the  impact  of  precipitation  particles  or  other  charging  mechanisms. 

The  discharges  create  avalanche  currents  in  the  air  and  launch  transient 
current  pulses  on  the  airframe  which  travel  from  the  discharge  point  to  the 
open  ends  of  the  airframe  where  they  are  reflected.  The  airframe  thus  be¬ 
haves  like  a  pulsed  transmission  line  and  radiates  whenever  the  current 


pulses  encounter  a  change  in  impedance  due  to  changes  in  the  geometrical 
or  electrical  characteristics  of  the  airframe. 


The  local  fields  produced  by  the  traveling  current  pulses  couple 
strongly  into  on-board  antennas  and  produce  p-static.  The  radiated  field 
components  are  negligible  by  comparison  and  are  not  usually  an  important 


source  of  p-static  receiver  noise. 

The  investigation  of  e-m  radiation  from  corona  discharges  has 


covered  a  broad  range  of  experimental  and  analytical  work.  The  work  has 
included  studies  of  the  radiation  from  corona  currents  in  air  and  from 
point  discharges  at  one  end  of  a  charged  conducting  cylinder,  field  measure¬ 
ments  of  the  natural  noise  environment  at  sea,  and  some  observations  of  sig¬ 
nals  associated  with  low  flying  commercial  aircraft.  ~ 
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I  I  , _ EXPERIMENTAL  WORK 

A.  R-f  energy  from  a  corona  point 

An  investigation  of  the  power  radiated  in  the  30  MHz  to  9  GHz 
frequency  range  by  corona  point  discharges  in  air  was  carried  out  at  the 
Cherry  Creek  Field  Site  Southeast  of  Denver  after  a  receiving  system  was 
assembled  from  available  components.  The  system  utilized  an  HP8555  Spec¬ 
trum  Analyzer  as  the  receiver  and  various  electric  field  probes  and  micro- 
wave  horn  antennas  as  sensors.  The  system  sensitivity  was  approximately 
-90  dbm  for  a  300  kHz  bandwidth,  the  widest  bandwidth  available  on  the 
HP8555  analyzer.  The  experiments  were  performed  inside  a  shielded  metal 
building  to  reduce  the  interference  from  commercial  radio  transmissions. 

R-f  energy  at  frequencies  below  150  MHz  was  found  to  be  detectable 
at  a  distance  of  1  meter  from  the  corona  point  with  this  recording  system. 
However,  no  signal  power  could  be  detected  at  higher  frequencies  due  to  the 
relatively  high  level  of  internal  system  noise  above  150  MHz. 

The  microwave  measurements  were  repeated  at  Fort  Huachucha, 
Arizona,  using  equipment  with  a  better  signal  to  noise  ratio  and  a  wider 
bandwidth  which  was  borrowed  from  ASATEC.  This  receiving  system  utilized  a 
Microtel  receiver  with  a  20  MHz  bandwidth  for  the  1  to  10  GHz  frequency 
range  and  a  Watki ns-Johnson  receiver  with  a  3  MHz  bandwidth  to  cover  fre¬ 
quencies  below  1  GHz.  Various  TWT  preamplifiers  were  used  to  Increase  the 
system  gain. 

R-f  energy  from  the  corona  point  was  found  to  be  detectable  with 
this  receiving  system  over  the  entire  frequency  range  studied.  The  results 
are  shown  in  Figure  1  as  a  function  of  frequency  with  the  measurements  below 
1  GHz  scaled  as  -/20/3  to  account  for  the  different  bandwidths  used.  Actual 
measurements  at  different  bandwidths  indicated  that  the  signal  amplitude  In¬ 
creased  with  bandwidth  somewhat  more  than  proportional  to  the  square  root  of 


bandwi dth. 


The  corona  current  was  monitored  during  the  experiments  by 
measuring  the  voltage  pulse  across  a  50  ohm  high  frequency  disk  resistor. 

A  typical  current  pulse  recorded  with  a  Tektronix  R7912  Transient  Digitizer 
with  a  rise  time  of  approximately  1  nanosecond  is  shown  in  Figure  2  For 
negative  point  corona. 

B.  Corona  from  the  end  of  a  cylindrical  conductor 

1.  Description  of  experiments 

Experiments  were  performed  on  10  and  20  foot  lengths  of 
2-1 A  inch  diameter  aluminum  pipe  with  the  following  objectives: 

a.  Verify  the  pulsed  transmission  line  model  of  the 
current  distribution. 

b.  Measure  the  corona  repetition  rate  and  d.c.  charging 
current  as  a  function  of  the  charging  voltage. 

c.  Measure  the  radiated  power  in  the  far  zone  (d  >  1.2  x 
wavelength)  of  the  cylinder. 

The  work  was  done  at  the  Cherry  Creek  Field  Site  with  the 
pipe  supported  horizontally  or  vertically  by  an  insulated  stand.  One  end 
of  the  pipe  was  capped  with  a  spherical  conductor  and  the  other  was  tapered 
to  a  point.  Voltages  of  up  to  -50  KV  were  applied  through  a  charging  coil 
so  that  any  transient  current  pulses  initiated  by  a  corona  discharge  at  one 
end  would  travel  cn  the  pipe  rather  than  through  the  power  supply. 

2.  Current  distribution 

The  current  distribution  on  the  pipe  was  measured  by  mount¬ 
ing  the  20  foot  pipe  horizontally  1  meter  above  a  wire  mesh  ground  screen 
and  recording  the  radial  electric  field  at  the  ground  directly  below  the 
pipe.  For  this  geometry  the  radial  electric  field  and  current  are  related 
by  the  equation ^ 


120  Kt) 
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where  h  Is  the  distance  to  the  ground  plane. 

The  E-field  measurements  were  made  with  a  short  probe  an¬ 
tenna  which  fed  a  high  Impedance  preamplifier,  followed  by  a  40  MHz  high 
pass  filter  to  reduce  the  background  noise  from  commercial  radio  transmis¬ 
sions.  The  output  was  recorded  on  film  with  a  Tektronix  454  oscilloscope. 
The  overall  system  passband  was  40  to  150  MHz  for  this  experiment. 

Measurements  were  obtained  with  the  probe  at  distances  of 
4,  8,  and  12  feet  from  the  corona  point.  The  results  are  shown  In  Figure  3, 
along  with  the  system  response  to  a  step  function  of  voltage.  The  initial 
pulse  and  reflections  from  the  open  ends  of  the  pipe  are  seen  to  occur  as 
expected.  The  reflection  time  delays  correspond  to  the  pulse  travel  time 
from  the  observation  point  to  the  ends  of  the  pipe  and  back  at  approximately 
the  speed  of  light.  The  system  sensitivity  was  such  that  a  field  strength 
of  3*3  volts/meter  corresponded  to  1  volt  deflection  on  the  oscilloscope. 

The 
5  m 
the 
the 

where  the  tin*  separation  is  maximum.  The  ratio  of  the  reflected  to  Initial 


;  170  mv  amplitude  of  the  first  pulse  corresponds  to  a  current  pulse  of 
III  lamps.  The  first  reflection  from  the  end  of  the  pipe  arrived  before 
initial  pulse  had  died  out  completely.  This  distorted  the  amplitude  of 

I 

I  reflected  pulse;  however  the  distortion  is  low  at  the  4  foot  position 


pulse  amplitude  at  this  point  gives  a  reflection  coefficient  of  approximately 
0.8.  The  reflection  coefficient  would  be  unity  in  the  Ideal  case. 

3.  Corona  pulse  rate  and  charging  current 

The  corona  pulse  rate  and  d.c.  charging  current  on  the  20  foot 
pipe  were  measured  at  pipe  voltages  between  -6  KV  and  -30  KV.  Trigger  pulses 
were  picked  up  from  a  small  metal  plate  placed  approximately  10  cm  from  the 
corona  point  and  grounded  through  a  50  ohm  resistor.  The  voltage  pulse 
across  the  resistor  caused  by  induced  current  in  the  plate  was  amplified  and 
used  to  trigger  an  electronic  counter. 


«>  CURRENT  DISTRIBUTION  12* 
FROM  THE  CORONA  POINTS 
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FIGURE  ^  CURRENT  DISTRIBUTION  ON  20'  PIPE 
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The  average  pulse  rate  as  a  function  of  the  charging  volt¬ 
age  Is  shown  In  Figure  4.  The  rate  increases  sharply  once  the  corona 
threshold  voltage  of  approximately  -5  KV  Is  exceeded. 

The  charging  current  and  average  current  per  pulse  as  a 
function  of  voltage  are  shown  In  Figure  5.  The  increase  in  the  average 
current  per  pulse  as  the  voltage  increases  is  consistent  with  an  expected 
increase  in  the  number  of  electron  avalanches  per  pulse  as  the  field  energy 
increases. 

4.  r-f  radiation 

The  power  radiated  by  vertical  10  foot  and  20  foot  lengths 
of  pipe  was  measured  at  distances  of  25  and  50  feet.  At  these  distances 
the  far  field  is  a  factor  of  10  greater  than  the  near  field  at  frequencies 
above  47  MH?  at  25  feet  and  24  MHz  at  50  feet.  The  measurements  were  made 
with  a  log-periodic  antenna  with  a  gain  of  12  db,  designed  for  the  30  to 
300  MHz  frequency  range.  The  Instrumentation  consisted  of  a  low  noise 
preamplifier,  150  feet  of  RG-8  coax  cable,  and  an  HP8555A  spectrum  analyzer. 
The  data  were  recorded  by  photographing  the  display  scope  of  the  spectrum 
analyzer  when  It  was  operated  In  a  swept-frequency  mode,  or  by  photographing 
the  detected  I.F.  output  on  an  auxilliary  scope  when  It  was  operated  as  a 
fixed  tuned  receiver.  The  system,  excluding  the  antenna,  was  calibrated  wtth 
a  VHF  signal  generator  over  the  10  to  480  MHz  frequency  range. 

Swept  frequency  recordings  of  the  power  received  In  a  300  kHz 
pass  band  from  the  10  foot  pipe  are  shown  in  Figure  6.  The  upper  trace  (a) 
shows  the  normal  background  noise  level  at  the  Cherry  Creek  Field  Site  with 
a  fairly  quiet  zone  between  approximately  25  and  60  MHz.  Trace  (b)  was 
recorded  at  25  foot  with  a  single  corona  point  at  the  upper  end  of  the  .  Ipe. 
The  lower  spherically  capped  end  was  24"  above  ground  level.  Traces  (c) 
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AVERAGE  CURRENT  PER  PULSE  -  NANOAMPS 


(«)  BACKGROUND  NOISE  («km) 


(I)  CORONA  NOISE  AT  23* 
SINGLE  POINT  AT  TOP 


(c)  CORONA  NOISE  AT  23* 
SINGLE  POINT  AT  BOTTOM 


(«)  CORONA  NOISE  AT  SO* 
SINGLE  POINT  AT  BOTTOM 
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FIGURE  ^  R-f  POWER  RECEIVED  FROM  A  10*  PIPE  IN  CORONA 
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and  (d)  show  the  noise  power  at  25  feet  and  50  feet  respectively  with  the 
corona  point  at  the  lower  end  of  the  pipe. 

The  pipe  potential  was  held  at  **25  KV  for  these  experiments. 
A  broad  spectrum  of  corona  noise  Is  seen  at  levels  of  15  to  20  db  above 
the  background  level  over  the  25  to  60  MHz  frequency  range  where  the  back¬ 
ground  is  low  (-75  dbm) .  The  reduction  in  peak  amplitude  of  approximately 
6  db  between  the  25  foot  and  50  foot  recording  distances  is  consistent  with 
the  expected  6  db  change  that  should  occur  in  the  radiation  zone  when  the 
distance  is  doubled.  No  significant  differences  are  observed  between 
traces  (b)  and  (c)  which  correspond  to  corona  from  the  top  and  the  bottom 
of  the  pipe. 

The  power  measurements  obtained  with  the  spectrum  analyzer 
are  related  to  the  incident  field  by  the  equation 


where  Ej  ■  spectral  amplitude  of  th.  incident  field  (v-sec/m), 

f  «  receiver  center  frequency  (Hz), 
o 

P  »  recorded  power  (Watts) , 
g  -  voltage  gain  of  the  receiving  system. 

Equation  (2-2)  is  based  on  the  assumption  that  the  spectral 
amplitude  Ej  is  constant  over  the  receiver  pass  band  and  that  the  antenna 
gain  was  12  db.  As  an  example,  a  recorded  power  level  of  -60  dbm  at  a 
frequency  of  50  MHz  corresponds  to  an  incident  field  of  3  x  10  v-sec/m*. 

Strong  interference  at  frequencies  above  60  MHz  prevented 
measurement  of  the  corona  spectrum  at  higher  frequencies  using  the  swept 
frequency  mode.  However,  with  the  receiver  tuned  manually  to  frequency 

*  The  system  gain  at  50  MHz  was  a  factor  of 
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points  were  the  inte  ference  was  relatively  low,  the  corona  noise  could 
be  detected  by  observing  the  detected  I.F.  output  of  the  spectrum  analyzer. 
With  this  technique,  corona  noise  could  be  observed  at  frequencies  up  to 
approximately  120  MHz.  Examples  of  the  detected  output  In  a  300  kHz  pass 
band  are  shown  in  Figure  7  for  center  frequencies  of  30,  50,  75,  and  100 
MHz.  These  data  were  recorded  at  a  distance  of  25  feet  from  the  20  foot 
pipe,  .oh  was  held  at  a  potential  of  -8  KV. 

The  data  presented  here  were  all  recorded  with  the  radiating 
structure  and  receiving  antenna  oriented  vertically  over  soil  of  average 
conductivity.  The  fields  are  therefore  approximately  doubled  over  those 
that  would  be  observed  from  the  same  structure  in  the  absence  of  the  ground. 
Some  attempts  were  made  to  record  the  corona  noise  with  the  antennas  parallel 
to  the  ground;  however,  the  ground  reflected  signal  in  this  case  nearly  can¬ 
cels  the  direct  signal  after  a  time  delay  which  is  proportional  to  the  ray 
path  difference  between  the  direct  and  ground  reflected  field  components. 

The  time  delays  for  the  source-receiver  geometries  that  could  be  practically 
attained  was  just  a  few  nanoseconds  and  no  data  were  recorded  from  the  hori¬ 
zontal  configuration.  This  problem  is  discussed  in  more  detail  in  Section 
III  of  this  report. 

C.  Radio  Noise  Survey 

l.  Description  of  experiments 

Experiments  designed  to  measure  the  natural  noise  background 
in  the  30  to  350  MHz  frequency  range  were  carried  out  aboard  the  R/V  Acania, 
a  research  vessel  operated  by  the  Department  of  Oceanography  of  the  Naval 
Postgraduate  School,  Monterey,  California.  The  primary  goal  was  to  record 
the  broadband  impulsive  background  noise  at  sea  with  a  secondary  objective 
of  measuring  background  noise  at  narrower  bandwidths  on  a  swept  frequency 
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The  broadband  measurement  system  consisted  of  a  short  verti¬ 
cal  monopole  antenna  which  fed  a  high  impedance  preamplifier,  pulse  ampli¬ 
fiers,  and  a  Tektronix  485  oscilloscope.  The  system  pass  band  was  30  MHz 
to  450  MHz.  The  narrow  band  swept  frequency  system  utilized  a  log-periodic 
antenna  and  a  HP8555A  spectrum  analyzer  which  covered  the  50  to  500  MHz 
frequency  range. 

2.  Results 

The  experiments  were  performed  in  July  1975  at  a  maximum 
distance  from  shore  of  120  miles.  Figure  8  shows  a  typical  swept  frequency 
plot  of  the  noise  received  in  a  300  kHz  pass  band  using  the  narrow  band 
system.  Man-made  noise  at  levels  of  -100  to  -120  dbm  predominated  at  this 
distance  from  shore  over  much  of  the  frequency  range  studied.  At  closer 
distances  the  noise  was  higher  and  consisted  solely  of  man  made  signals. 

Measurements  of  the  impulsive  background  noise  were  unsuccess¬ 
ful  because  of  the  limited  distance  from  shore  of  the  Acania.  Television 
and/or  FM  broadcast  station  signals  at  the  maximum  distance  of  120  miles 
were  still  much  stronger  than  the  impulsive  background  noise  at  this  distance. 
D.  Field  Measurements  on  Low  Flying  Aircraft 

1.  Description  of  Experiments 

The  measurement  of  noise  generated  by  aircraft  in  flight 
was  attempted  with  available  instrumentation  after  modifying  the  HP8555A 
spectrum  analyzer  to  give  a  10  MHz  i.F.  bandwidth.  The  recording  system, 
using  the  log-periodic  antenna,  was  set  up  at  a  distance  of  100  meters  South 
of  the  East  end  of  the  East-West  runway  at  Stapleton  International  Airport 
in  Denver.  Observations  were  made  on  several  occasions  during  varied 
weather  conditions.  The  analyzer  was  tuned  to  relatively  quiet  frequency 
spots  near  150  MHz  or  275  MHz  and  the  i.f.  output  was  monitored  on  an  os¬ 
cilloscope. 
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2.  Results 

Signals  consisting  of  bursts  of  random  appearing  pulses, 
often  accompanied  by  a  short  quasi-CW  burst  were  observed  on  several 
occasions  as  aircraft  flew  through  the  antenna  beam.  The  most  consistent 
observations  were  obtained  during  a  two  hour  period  when  light  rain  was 
failing.  Nearly  all  of  the  aircraft  gave  detectable  signals  during  this 
period,  whereas  only  the  larger  aircraft  such  as  DC-10's  or  7^7's  were 
detectable  during  balmy  weather.  It  was  observed  that  automobile  ignition 
noise  was  generally  detectable  at  considerably  greater  distances  than  were 
aircraft. 

Although  these  experiments  Involved  a  large  number  of 
variables  whidh  could  not  be  controlled,  strong  correlation  was  observed 

between  the  passage  of  an  aircraft  through  the  antenna  beam  and  the  noise 

| 

output  of  the  recording  system.  However,  the  recorded  noise  could  not  be 

I 

definitely  attributed  to  corona  as  no  independent  method  was  available  to 

determine  whether  or  not  the  aircraft  were  actually  In  corona  at  the  time 

i 

the  observations  were  made. 
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III. 


ANALYSIS 


A.  The  corona  current  in  air 

A  corona  discharge  from  a  negatively  charged  conductor  produces 
electronic  and  Ionic  currents  in  the  air  near  the  corona  point  by  electron 
avalanching.  The  current  Is  dominated  Initially  by  the  motion  of  the  elec¬ 
trons  and  tatter  by  the  build  up  of  heavy  positive  and  negative  Ions. 
Photons  emitted  by  the  corona  envelope  cause  further  Ionization  of  the  atr 
and  also  release  new  electrons  from  the  conductor  by  photoemlssion.  The 
net  result  Is  a  very  fast  rising  pulse  of  current  during  the  electronic 
phase  which  decays  exponentially  In  time  during  the  Ionic  phase  of  the  dis¬ 
charge.  Only  the  electronic  current  component  has  been  studied  here,  since 
it  is  much  faster  and  contributes  most  of  the  high  frequency  content  of  the 
radiated  pulse. 

The  electron  current  density  in  air  was  calculated  by  solving 

Townsends  equations  for  an  electron  In  a  non-uniform  electric  field  between 

(2) 

a  negatively  charged  point  and  a  grounded  plane. 

The  result  Is  given  by  the  equation 

s  qryjS(  t-r)axp{(Ap/C,  )[l-exp(-C,X')]axp(-Cla)}  3.1 

where 

•  tq 

q  ■  electron  charge  (1.602  x  10  ^  coulombs) 

2 

nQ  ■  number  density  of  Initial  electrons  (electrons/an  ) 

A  ■  constant  (l1*. 6) 
p  ■  atmospheric  pressure  (760  mm) 

C|  •  constant  (160) 
a  "  rad -us  of  corona  point  (0.01  cm) 


The  parameter  t  was  determined  from  the  differential  equation  for  the  elec¬ 
tron  velocity  and  from  the  mobility  equation,  which  relates  the  velocity  to 


the  electric  field.  The  solution  gives 

-  Ci<2ox*-k|2) 

T=  2Bpb 


where 


3.2 


b  «  electron  mobility  (500  cm/sec  per  volt/cm), 

B  -  constant  (365). 

B.  The  corona  current  on  a  conducting  cylinder 

The  pulse  of  current  launched  onto  a  conductor  when  a  corona  dis¬ 
charge  occurs  can  be  determined  by  equating  the  work  done  by  the  field  to 
the  change  In  energy  of  the  charged  conductor.  This  gives 

Kt'j^/EU’N.U'.t'jdV  3.3 

where  l(t')  Is  the  current  pulse  on  the  conductor.  This  equation  was  evalu¬ 
ated  by  substituting  equation  3*1  for  Jj  and  approximating  the  field  distri¬ 
bution  by  the  equation 


E(*'l 


3.4 


(a  +  x')ln(2i/a) 

which  is  a  sufficient  approximation  for  the  field  between  a  point  of  radius 
"a"  raised  to  a  potential  V  and  a  plane  at  a  distance  t  from  the  point. ^ 
With  these  substitutions  the  Integral  in  equation  3*3  evaluates  to 
_  |0a«p((Ap/C,)[l-axp(“C,VQ^2Bpbtl/C,-o)axp(-Ctq)]} 


MO 


a*+  2Bp'bt’/C, 


where 


l8.  2Ld*q"0BP,b 

C,  In  (2 |/a) 


3.5 


3.6 


The  corona  channel  was  assumed  to  be  cylindrical  of  radius  d(cm)  In  this 
derivation.  The  current  pulse  given  by  equation  3.5  la  shown  plotted  In 
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Figure  9.  using  appropriate  values  for  the  various  constants  In  equations 

(M 

3.5  and  3.6*  consistent  with  a  typical  discharge. 

C.  The  radiated  fields 


The  fields  radiated  by  the  corona  current  In  air,  and  the  current 


on  the  cylinder  have  been  calculated  separately  by  evaluating  the  vector 
potential  for  each  case.  The  vector  potential  is  given  by  the  equation ^ 


V'.*  1  * 


where  the  subscript  I  applies  to  the  corona  current  In  air  and  the  subscript 
2  applies  to  the  current  on  the  cylinder.  The  constant  c  Is  the  speed  of 
light  and  the  primed  and  unprimed  variables  -efer  to  the  source  and  receiver 
coordinates  respectively.  The  distance  r  Is  measured  from  the, current 
element  J(x',t')  to  the  receiving  point.  The  geometry  Is  shown  In  Figure  10. 

i 


Figure  10.  Source  Geometry 


The  electric  field  component  E0  In  the  radiation  zone  Is  related 
to  the  vector  potential  by  the  equation 

E-«.i  3.8 

^  e  at  * 


1  V  WW  VJiV 


•  Vt.  >  v « V*  •*>  v»  *.*•  *>  v*  »>  ; 


TOTAL  CURRENT 


FIGURE  9.  CALCULATED  CURRENT  PULSE  ON  A  CYLINOER 
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and  the  magnetic  field  H.  Is  then  given  by  H.  «  E./{120ir).  3.9 

t  f  o 

The  integration  of  equation  3.7  was  carried  out  directly  In  the 
time  domain  In  region  1,  using  equation  3.1  for  Jj.  The  final  result  can 
be  written  in  the  form 


where 


m  *  2Bpb/C, 


3.10 


3.1! 


C*s-Ap/C,  3.12 

Cs=  vrd*qm  exp  [-C2exp(- C,o)]/2  3.13 

and  the  approximations 

r  ssr  R  ( l-x'coe Q/R)  3.1* 

dV«Vdadx*  3.15 

x/R«l  3.16 


have  been  used.  Equation  3.10,  which  gives  the  electric  field  pulse 
radiated  by  the  corona  current  In  air  is  shown  plotted  In  Figure  11. 

The  fields  radiated  by  the  current  on  the  cylinder  In  region  2 
were  computed  In  the  frequency  domain  and  then  transformed  beck  Into  the 


time  domain.  For  this  case  the  vector  potential  has  the  form 

...  .  KtwIexpMkRj/dx'Kx'jexpHkx'eostf) 

A(R,w)  - - fr - - 


FIGURE  II.  RADIATED  PULSE  FROM  CORONA  ELECTRON  CURRENT  IN  AIR 
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where  K(u)  is  the  spectrum  of  the  current  pulse  defined  by  equation  3-5 
and  plotted  in  the  time  domain  in  Figure  9* 

I(x')  Is  the  current  distribution  along  the  cylinder  and  k  ■  u/c 
is  the  propagation  constant  in  free  space.  I(x')  can  be  determined  from 
transmission  line  theory,  which  for  a  thin  cylinder  has  the  form^ 

I  (x‘)  =  exp(-yx') +  T,  axp[-y(2L-x')]  +  T0T,  exp[-y(2L+w  }]  3.18 

♦  T0T,2 exp  [-y (4 L-x '  T| T,2  exp  [- y  ( 4 L x' )] 

♦TjT’exp  [-y(6L-x‘  )]+•- 

The  substitution  of  equation  3.18  into  3.17  and  evaluation  of  the  integral  • 
gives  the  vector  potential  In  region  2.  The  radiated  electric  field  «»t 
point  P  is  then  determined  from  equation  3.8.  The  result  Is 

E(Rlaj,»^^.exp(-ikR)tiw)K(w,{~~L  [exp  Hkz,) 

♦T0T,  exp  (-IM2U- *,))+•  •  ]  3*!9 

♦*^[T,  exp(ik  (2L-z2  ))  +  T0T*exp  I- ik(  4L- z2  ))♦•  •  •]} 


The  above  results  pertain  to  a  cylinder  of  length  L  where  Tq  Is  the  reflec¬ 
tion  coefficient  at  the  corona  end  of  the  cylinder  and  T|  Is  the  reflection 
coefficient  at  the  opposite  (open)  end.  y  is  the  transmission  constant 
given  by 

y -  a  +ik 


/ 


3.20 
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where  e  gives  the  attenuation  and  k  is  the  phase  constant  of  the  cylinder 
when  treated  as  a  transmission  line.  The  variables  Zj  and  Z^  are  defined 
by  the  equation 

Zl|2*  L{  I  ±  cos  9)/2  3.21 

For  the  special  case  8  -  n/2,  Zj  ■  Zj  ■  L/2  and  equation  3«19  can  be  written 
In  MKS  units  as 

E(R,bi)  =  '-(sin  irL/\){exp  [’ik(R  +  L/z]  +  •  •  •}  3.22 

where  the  terms  In  braces  correspond  to  the  successive  reflections  from  the 
two  open  ends  of  the  cylinder. 

The  inverse  transform  of  equation  3 - T9  gives  the  time  domain 
solution.  This  Is  most  easily  obtained  by  using  the  time  convolution 
theorem  for  the  products  of  frequency  functions,  which  has  the  notation 

•  •  —  F,  (u»)*  Fj  (o») ■  •  •  •  3.23 

where  the  asterisk  denotes  convolution  and  the  arrow  signifies  the  Fourier 
transform.  The  final  result  In  MKS  units  can  be  written 


tf(R,t)«  -^-i(D*{^j[8(t*)-8(t#-2Z|/C)+T0Tl8(t#-2LA:)3.2l, 

•  T#T,  8(»*  2L/C-2Z,/C  )  ♦  •  •  •] 

+  [t,  8  (t  ■  2L/C  +  2Zt/C)-T,8(t*  2L/CK-.  •]} 

where  t*  ■  t  -  R/C. 


'  /\ 

;  /  \\V' 
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The  reflection  coefficients  Tg  and  Tj  have  been  assumed  Indepen¬ 
dent  of  frequency  In  equation  3*24,  so  that  only  the  amplitude  and  not  the 
shape  of  the  current  pulse  I (t)  Is  changed  upon  reflection. 

The  form  of  equation  3.24  Is  convenient  for  computations  because 
the  convolution  of  the  current  pulse  I (t)  with  the  string  of  delta  functions 
contained  within  the  braces  simply  superimposes  the  pulses,  each  shifted  in 
time  to  that  given  by  the  argument  of  each  delta  function.  The  results  for 
the  special  case  6  ■  ir/2  and  Tg  ■  T^  ■  -0.8  are  depicted  dlagramatlcal ty  In 
Figure  12. 

All  of  the  calculations  have  been  carried  out  for  a  conductor  In 
free  (air)  space.  The  results  for  a  conductor  near  the  earth  would  be  modi¬ 
fied  by  the  Image  currents  In  the  ground.  The  effect  on  the  radiated  signal 
from  a  vertical  conductor  over  ground  of  finite  conductivity  would  be  to  In- 
.crease  the  signal  amplitude  by  the  factor  (HR),  where  R  Is  the  ground  re¬ 
flection  coefficient,  which  Is  positive  for  vertical  polarization.  For  In¬ 
finite  conductivity,  R  •  +1  and  the  signal  amplitude  would  be  doubled. 

The  horizontally  polarized  field  radiated  by  an  elevated  conduc¬ 
tor  mounted  parallel  to  the  earth's  surface. can  also  be  estimated.  For 
this  case  the  ground  reflection  coefficient  Is  approximately  equal  to  -1 
for  all  frequencies.  The  resultant  field  at  an  elevated  receiving  antenna, 
which  Is  the  sum  of  a  direct  and  a  ground  reflected  component,  would  be 
proportional  to  the  derivative  of  the  free  space  field  over  a  restricted 
frequency  range.  This  can  be  seen  from  the  equation 

E(t)  «  *(»)«-R*(t-r)  3.25 

where 

e(t)  Is  the  free  space  field, 

R  Is  the  reflection  coefficient  for  horizontal  polarfzatlon(R*  -1), 

T  Is  the  path  difference  delay  time. 

E(t)  Is  the  resultant  field  at  the  antenna. 
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The  substitution  t'  ■  t  +  t/2  and  R  ■  -1  gives 
E(t‘)  =  Q  (t'+r/2)-tf  ( t‘-r/2)  3.26 

which  has  the  Fourier  transform 

6(<u)  =  loir  F(tu)  3.27 

where  G(u)  Is  the  Fourier  transform  of  E(t)  and  F(u)  Is  the  Fourier  trans¬ 
form  of  e(t).  When  wx  Is  small, 

6(ui)ofl(urF(tu)  3.28 

which  transforms  back  Into  the  time  domain  as 

E{f)stTMi!  3.29 

d  i 

This  shows  that  the  resultant  horizontal  field  at  the  antenna  Is  approxi¬ 
mately  proportional  to  the  derivative  of  the  free  space  field.  The  fre¬ 
quency  range  over  which  3*27  Is  valid  Is  determined  by  the  sine  function. 

Thus  for  10%  accuracy  the  requirement  Is  • 

0.9  <  sin  c(orr/2)  S  I  3.30 

<  • 

which  Is  met  when  fir  -  0.25.  As  an  example,  If  the  time  difference  x  was 

5  nanoseconds,  the  resultant  field  would  be  proportional  to  the  derivative 
(within  10%) for  frequencies 

f  -  0. 25/(5  x  lo"9), 
or 

f  -  50  MHz. 

An  order  of  magnitude  estimate  of  the  radiated  spectral  amplitude 
E(R,u)  can  be  obtained  from  equation  3>22  once  the  spectrum  K(u)  of  the 
current  pulse  is  known.  This  could  be  computed  numerically  from  the  recorded 
corona  current  pulse  shown  In  Figure  2;  however,  for  the  frequency  range  of 
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Interest,  the  pulse  Is  adequately  represented  In  closed  form  by  the  double 
exponential  function 

l(t)  *  l0  [axp(-ot)  -exp(-bf)]  3.31 

which  has  the  Fourier  transform 

-  r.  J?r-T»  352 

where  I  Is  the  peak  amplitude  of  the  current  pulse  and  the  decay  constants 
o 

a  and  b  have  the  values 

a  -  3.3  x  107,  b  ■  9.4  x  109.  With  l#  ■  5  «,  f  ■  50  MHz, 
the  magnitude  of  the  spectrum  Is 

| K ( f  )|  *  1.6x10  '*  amp-sac  3*33 

Substitution  of  this  value  Into  equation  3.22  gives 

E(f)  »  2  x  10“9  v-sec/m  at  I  meter  for  a  single  pulse  train 
radiated  by  a  cylinder  of  length  L  •  10  feet. 


*** 


•  N 

r 


IV. 


CONCLUDING  REMARKS 


The  radiation  from  a  conducting  cylinder  in  corona  at  one  end 
has  been  measured  and  found  to  be  in  good  agreement  with  calculations  of 
the  current  pulse  which  travels  on  the  cylinder  and  r’  transmission  line 
model  of  the  current  distribution.  As  an  example,  the  -56  dbm  power  level 
measured  at  50  MHz,  with  the  receiving  antenna  at  a  distance  of  25  feet 

•  Q 

from  the  10  foot  pipe  corresponds  to  a  spectral  amplitude  of  3.6  x  10 
v-sec/m  at  1  meter.  This  compares  very  well  with  the  calculated  value  of 

.A 

approximately  4  x  10  J  v-sec/m  for  a  vertical  cylinder  over  a  perfect  ground 

plane.  For  a  random  discharge  rate  of  100  kHz,  the  incident  flux  density 

-15  2 

would  be  approximately  4  x  10  s  watts-sec/m  compared  to  a  background  noise 

“  1  Q 

flux  from  atmospheric  and  galactic  noise  of  approximately  1.4  x  10  J 

2  (7) 

watts-sec/m  in  the  50  to  100  MHz  frequency  range. 

Cylinders  were  used  in  these  experiments  because  the  simple  geo¬ 
metry  involved  is  amenable  to  calculation.  More  complicated  structures  such 
as  airframes  would  be  expected  to  behave  differently;  however,  the  basic 
transmission  line  model  is  directly  applicable  to  more  complicated  geometrical 
shapes.  Scaled  models  of  aircraft  could  be  used  to  determine  many  of  the 
characteristics  of  the  radiated  fields  expected  from  full  scale  aircraft  in 
corona. 

The  charging  current  and  potential  of  the  cylinders  studied  were 

(0) 

much  lower  than  the  nominal  values  given  for  aircraft  In  flight.  '  It  is 
reasonable  to  expect  that  the  incident  noise  flux  from  an  airplane  in  corona 
would  be  considerably  more  Intense  than  the  flux  which  has  been  determined 
from  these  experiments  with  cylinders. 
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